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(BESIII Collaboration) 

The number of J ftp events collected with the BESIII detector at the BEPCII from June 12 to 
July 28, 2009 is determined to be (225.3 ± 2.8) x 10® using J/ip inclusive events, where the 
uncertainty is the systematic error and the statistical one is negligible. 



INTRODUCTION Electron- Positron Collider (BEPC) and the Beijing Spec- 



To meet the challenge of precision measurements 
of T— charm physics, a major upgrade on the Beijing 
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trometer (BES) was completed in 2008 (now called 
BEPCII and BESIII). BEPCII is a double ring e+e" col- 
lider with a design peak luminosity of 10^'^ cm~^s~^ at 
^/s =3.773 GeV, which is 100 times that of its predeces- 
sor. The BESIII detector is a large solid-angle magnetic 
spectrometer that is described in detail in Ref. The 
major improvements in the BES detector are the huge 
superconducting solenoid magnet with a central field of 
1 T, which offers a significant improvement in the mo- 
mentum resolution of charged particles, and a cesium 
iodide (Csl) calorimeter for the energy measurement of 
electrons and photons, which provides more than a fac- 
tor of 10 improvement in the precision of electromagnetic 
shower energy measurements. 

Since the discovery of the J/V' in 1974, it has al- 
ways been regarded as an ideal laboratory to study 
light hadron spectroscopy and to search for new types 
of hadrons (e.g. glueballs, hybrids and exotics). With 
58 million J/ if} events collected with the BESII detec- 
tor, many important results have been obtained, which 
underlined the importance of the study of J /ip decays. 
Therefore, after a successful commissioning of the BE- 
SIII detector together with BEPCII, a large sample of 
J /ip events was collected from June 12 to July 28, 2009, 
which allows the study of the properties and the decays 
of the J/V' with unprecedented precision. 

The number of J /tp events and its uncertainty are two 
key quantities in the precision measurements of J de- 
cays. At BESII, the number of J/ip events was deter- 
mined with J /tp — >■ 4-prong events, and its systematic 
uncertainty was 4.7% The excellent BESIII detec- 
tor and its good performance allow the determination of 
the number of J /ip events with higher precision. To re- 
duce the systematic uncertainty from that in Ref. 0, a 
new method using J /ip ^ inclusive events is introduced. 
The number of J/'ip events {Nj/^) is calculated with 



j/i, 



(trig X e^^(^ X /, 



(1) 



where Nsei is the number of J/ip — > inclusive events se- 
lected from J/ tp data; Nbg is the number of background 
events estimated from the continuum data taken at the 
center-of-mass energy of 3.08 GeV; etrig is the trigger 

efficiency; e^^^^ is the J/ip — >■ inclusive detection ef- 
ficiency determined experimentally from ip' data using 
ip' — > 7r+7r~J/'0 events; fcor is a correction factor for 
^data' obtained from Monte Carlo (MC) simulation which 
accounts for the difference between the J/tp events pro- 
duced at rest and those produced from tp' tt+tt" J/-0. 
The correction factor in Eq. ([l}, which is approximately 
unity, is determined from 



fcor — 



(2) 



where eilc is the detection efficiency of J/ip — 
inclusive events determined from the J /ip MC sample 



and ej^^ is the efficiency determined from the Tp' — 
TT^TT~ J/ip {J/ip — >■ inclusive) MC sample. 

There are two major improvements over the method 
in Ref. One is the generalization of the J/ip — 

4-prong events to J/ip — > inclusive events, which al- 
lows the number of J/ip events to be determined by re- 
quiring different numbers of charged tracks; the other 
is to use the MC samples of J/ip — ^ inclusive and 
Ip' — > TT^TT^J/ip (J/ip — > inclusive) events generated 
with the BesEvtCen generator (31 based on GEANT4 [i| 
to determine the correction factor, fcor- In this analysis, 
the events with more than one charged tracks are used 
to determine the number of J/ip events. 

At present only about 50% of the J/ip decays are 
observed and listed in the Particle Data Group tables 
(PDG) Q. In the MC simulation package, the unknown 
J/ip decays are roughly generated with the Lundcharm 
model. In the Lundcharm model, charmonium decay via 
gluons is described by the QCD partonic theory, and the 
partonic hadronization is handled by the LUND model. 
Extended C- and G-parity conservation are assumed and 
abnormal suppression effects of charmonium decay are 
included [61] . 



II. 



J/iP 



inclusive SELECTION CRITERIA 



Event selection criteria are required to distinguish 
J/ip — inclusive events from Bhabha (e^e^ — !■ e+e^), 
dimuon (e"*"e~ — )■ j.i'^fi'^), cosmic ray and beam-gas 
events in J/ip data. 

At the track level, candidate events are required to 
satisfy the following selection criteria: 

1. Charged tracks are reconstructed using hits in the 
Main Drift Chamber (MDC) and are required to 
be in the polar angle range |cos0| < 0.93, have 
momentum p < 2.0 GeV/c, and have the point of 
closest approach of the track to the beamline within 
15 cm of the interaction point along the beam di- 
rection (Vz) and within 1 cm in the plane perpen- 
dicular to the beam (Vr). 

2. Clusters in the electromagnetic calorimeter (EMC) 
must have at least 25 (50) MeV of energy in the 
barrel (end cap) EMC, have | cos 6*1 < 0.83 in the 
barrel (0.86 < |cos0| < 0.93 in the endcap), and 
have EMC cluster timing T in the range of < T < 
15 (with unit of 50 ns) to suppress electronic noise 
and energy deposits unrelated to the event. 

At the event level, at least two charged tracks are 
required, and the visible energy, E^is, must be greater 
than 1.0 GeV. Here Ems is defined as the sum of charged 
particle energies computed from the track momenta by 
assuming pion masses and the neutral shower energies 
deposited in the EMC. According to the distribution of 
visible energy shown in Fig. [1] this requirement removes 
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two thirds of background events, estimated with the con- 
tinuum data taken at the center-of-mass energy of 3.08 
GeV, while it has httle effect on the inclusive events. 

To remove background from Bhabhas and dimuons, 
events with only two charged tracks must have the mo- 
menta of both charged tracks less than 1.5 GeV/c. Fig. [2] 
displays the scatter plot of the momenta of two charged 
tracks, where the clear cluster with the momenta around 
1.55 GeV/c corresponds to the contribution from lep- 
tonic pairs. Most of the leptonic pairs are removed by 
the above requirement as indicated by the solid lines in 
Fig.[5J From the deposited energy distribution of charged 
tracks in the EMC, shown in Fig. [3l a peak around 1.5 
GeV is clearly observed, which corresponds to the con- 
tribution of Bhabha events. Therefore, to further re- 
move Bhabha events, the deposited energy in the EMC 
of each charged track is required to be less than 1 GeV. 
After the momentum and energy selections there remain 
174.28 ± 0.01 million events (Nsei) from the J/ip data. 
The distributions of the track parameters for closest ap- 
proach and track angle Vr, Vz, cos 9, 

the total energy deposited in the EMC (Eemc), and 
the charged multiplicity (Ngood) after subtracting back- 
ground events estimated with the continuum data taken 
at the center-of-mass energy of 3.08 GeV (see Section Hill 
for details) are shown in Figs. U] through [51 respectively. 
Also shown are the distributions from MC simulation, 
normalized to J/ ip data. The distributions of Vz, Vr, and 
cos 9 of charged tracks, and the Eemc distribution for MC 
simulation are in reasonable agreement with those from 
data. 

For the charged multiplicity distribution shown in 
Fig. Ill neither the MC simulation with the Lundcharm 
model nor the MC simulation without the Lundcharm 
model agree very well with the data. However the effect 
of this discrepancy between data and MC simulation on 
the correction factor is very small, as described in Sec- 
tion EIIl 



III. 



BACKGROUND ANALYSIS 



Background events come mainly from Quantum 
Electro-Dynamics (QED) processes, beam-gas interac- 
tions, and cosmic rays. In this analysis, all of them are 
estimated with the number of events selected from the 
continuum data taken at the center-of-mass energy of 
3.08 GeV, normalized to the J/i/j data after taking into 
account the energy-dependent cross section of the QED 
process: 



bg 



N: 



£ 



3.08 X 



£ 



X . 



3.08 



(3) 



where Nbg is the estimated number of background events 
in the selected J/ip events; iV3.08 is the number of events 
selected from the continuum data; £j/^ and £3.08 are 
the integrated luminosities for J/ip and continuum data, 
respectively; ^sjj^ and ^sa.os are the center-of-mass 
energies for J jt^ data (3.097 GeV) and the continuum 
data (3.080 GeV). 

The integrated luminosities are determined using 
e~ — > 77 events with the following selection crite- 
ria: there are at least two neutral tracks with the de- 
posited energy of the second most energetic shower larger 
than 1.2 GeV and less than 1.6 GeV; and |cos6'| < 0.8, 
where 9 is the polar angle in the EMC. The number 
of signal events is determined by counting in the sig- 
nal region |A0| < 2.5° and the background estimated 
in the sideband region 2.5 < |A</i| < 5°, where A0 = 
(/).y2| — 180° and ^ is the angle of photon in x-y plane. 
Figs. Eland [TOlshow the distribution of energy deposited 
in EMC and cos 9 of photons. The integrated luminosi- 
ties of J/V' data and continuum data are determined to 
be 79631 ± 70 {stat) nb-^ and 281 ± 4 (stat.) nb'^, re- 
spectively. Here, the statistic error is 1.5%, and the sys- 
tematic error can be cancelled according to Eqs. [31 
With the same selection criteria for inclusive events 
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FIG. 1. The visible energy distributions for J/tp data (dots 
with error bars), continuum data (circles with error bars) and 
MC simulation of J/i/) — >■ inclusive events (histogram). The 
arrow indicates the minimum Evis required to select inclusive 
events. 



FIG. 2. The scatter plot of the momenta of the charged 
tracks for 2-prong events. The cluster around 1.55 GeV/c 
corresponds to the contribution from leptonic pairs. Most are 
removed with the requirements on the two charged tracks, 
Pi < 1.5 GeV/c and p2 < 1.5 GeV/c, as indicated by the 
solid lines. 
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FIG. 3. The distributions of deposited energy in the EMC by 
the charged tracks of 2-prong events for J/tp data (dots with 
error bars) and for the combined, normalized MC simulations 
of e'^e~ — > e'^e~ and J/ip — >■ e^e~ (histogram). 



from J/ip data, 21266 ± 146 events are selected from 
the continuum data. Therefore the number of back- 
ground events (Ntg) is estimated to be 5.96 ± 0.04 mil- 
lion using Eq. ([3]). The background ratio in the selected 
J/V' — >■ inclusive events is calculated to be 3.5% by com- 
paring the number of background events to the number 
of inclusive events selected from J/i/j data. 

In the above calculation, the background events from 
cosmic rays and beam-gas interaction are normalized 
with the same procedure as QED events. In fact, the 
number of cosmic rays is proportional to the data tak- 
ing time, whereas beam-gas events are related with the 
vacuum status and the beam current for taking data, in 
addition to the data taking time. In this analysis, the 
difference of the number of background events estimated 
with and without considering the energy dependence of 
the cross section for QED processes is taken into account 
in the overall systematic uncertainty of the number of 
J/ip events (see Section IVlII for details). 



IV. DETERMINATION OF THE DETECTION 
EFFICIENCY AND CORRECTION FACTOR 



Usually the detection efficiency is determined using a 
MC simulation of J ftp — > inclusive^ assuming that the 
detector response is well simulated. The efficiency is then 
the ratio between the number of events detected and the 
number of events generated. In this analysis to avoid the 
uncertainty caused by any discrepancy between MC sim- 
ulation and data, the detection efficiency is determined 
experimentally using 106 million -0' events taken with the 
BESIII detector. The experimental detection efficiency, 
^tata' then the number of selected events divided by 
all J/ip ^ inclusive events obtained from the cascade 
decays of ip' — tt+tt" J/-0. 



To select -0' 



TT J/ip events, there must be at least 



two soft pions that are each reconstructed successfully in 
the MDC within the polar angle range | cos 6*1 < 0.93, 
have Vr < 1 cm and \Vz\ < 15 cm, and have mo- 
mentum less than 0.4 GcV/c. The tt momentum dis- 
tributions in Fig. [11] show that the MC simulation is 
in good agreement with data. There are no other re- 
quirements on the remaining charged and neutral tracks. 
The invariant masses recoiling against all possible tt+tt^ 
pairs are calculated and shown in Fig. [T^l A clear 
peak around 3.1 GeV/c^, corresponding to the decay of 
Tp' — >■ 7r"'"7r~J/0, J ftp — ^ inclusive, is observed over a 
large flat background. The number of J/ip — inclusive 
events, N^c — (19526 ± 10) x 10^^, is obtained by a fit to 
the 7r"'"7r~ recoil mass spectrum with a double-Gaussian 
plus a second order Chebychev background function. 

To determine the number of selected J/ip — >■ inclusive 
events, in addition to the above common selection cri- 
teria for the two soft charged pions, the remaining 
charged tracks and neutral tracks must satisfy the re- 
quirements for the J/f/) — inclusive events described in 
Section IIIII Fig. [T^] shows the invariant mass recoiling 
against tt+tt^ for the selected events, and the number of 
selected J/0 inclusive events, Nf^j., is determined to 
be (14432±9) x 10'^ from a fit with a double-Gaussian plus 
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FIG. 4. The distributions of Vz for J/ip data (dots with error 
bars) and MC simulation of J/tp ^ inclusive (histogram). 



FIG. 5. The distributions of Vr for J/tp data (dots with error 
bars) and MC simulation of J/ip ^ inclusive (histogram). 
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a second order Chebychev background function. Finally 
the experimental detection efficiency of J /ip — >■ inclusive 

events, e^^ta^ determined to be (73.91 ± 0.02)%. 

Since the J/ij} decays in flight, a correction factor de- 
fined as in Eq. ([2]) is used to correct for the kinemati- 
cal effect in order to determine the detection efficiency 
for direct e^e~ — >■ J/ip — >■ inclusive decays. With 
the same procedure, including the event selection cri- 
teria and the fit functions, the detection efficiency of 
ej^'^ = (75.87 ± 0.06)%, is obtained from a MC sam- 
ple of 2 million of i/'' — ?> tt+tt^ J/-)/) events. To deter- 
mine eilc , a MC sample of 1 million events of J/ijj — > 
inclusive was generated. With the same selection crite- 
ria for J /ip — inclusive events as listed in Section m 
766893 ± 423 events are selected, and the corresponding 
detection efficiency is calculated to be (76.69 ± 0.04)%. 
The correction factor fcor for the detection efficiency, is 
then determined to be 

fcor = ^ = 1-0108 ± 0.0009. (4) 

V. TRIGGER EFFICIENCY 

The trigger efficiency of the BESIII detector has been 
studied using different physics channels [tI and was found 
to be very close to 100%. Therefore, we do not repeat a 
similar study here, but assume a 100% trigger efficiency. 

VI. THE NUMBER OF J/ip EVENTS 

The values of different parameters used in Eq. ([l} are 
summarized in Table HI and the number of J/ ip events 
is then calculated to be (225.30 ± 0.02) x 10^. Here the 
statistical error is only from Ng^i, while the statistical 
fluctuation of N^g is taken int account as part of the 



systematic uncertainties (see subsection 7.4). The sys- 
tematic errors from different sources will be discussed in 
the next section in detail. 



VII. SYSTEMATIC UNCERTAINTY 
A. MC model uncertainty 

The efficiency correction factor (fcor), which is used to 
correct the detection efficiency for the in-flight J/ ip decay 
from -0' data, is a MC simulation dependent parameter. 

To check the MC model dependence of the correction 
factor, we also determine the correction factor with MC 
samples generated without the Lundcharm model. The 
difference of the correction factors obtained with and 
without the Lundcharm model, 0.49%, is taken as the 
systematic uncertainty from the MC model in the deter- 
mination of the number of Jf-tjj events. 



B. Tracking efficiency 

According to tracking efficiency studies, the consis- 
tency of tracking efficiencies between MC simulation and 
data in J/tp decays is 1% for each charged track, al- 
though it is a little larger at low momentum. 

In this analysis, the consistency of tracking efficiency 
between MC simulation and data in -0' decays is assumed 
to be the same as that in J/^ decays. Actually there may 
be a difference in the two data sets taken at different 
center-of-mass energies. To estimate the corresponding 
uncertainty, the tracking efficiency in the J/ip MC sam- 
ple was varied by —0.5% for the tracks with momentum 
greater than 350 MeV/c and —1.0% for the tracks with 
momentum less than 350 MeV/c. The change of the cor- 
rection factor due to this variation leads to a change of 
0.40% in the number of J/i/j events, which is taken as the 
systematic uncertainty due to the tracking efficiency. 



25000 : 



-g 20000 ; 

H 15000 ; 



10000 r 
5000 r 
' 



-1.0 



COS 9 




(GeV) 



FIG. 6. The cos6' distributions of charged tracks for J/^p data 
(dots with error bars) and MC simulation oi J/ip inclusive 
(histogram). 



FIG. 7. The distributions of the total energy deposited in 
the EMC oi J/ip ^ inclusive events for J/t/; data (dots with 
error bars) and MC simulation (histogram). 
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C. Fitting of J/V' peak 

From the fit of tfie J/i/; peak we obtain tlie fitting er- 
rors 0.03% and 0.08% in tire determination of e^^jji and 
^mci respectively, In addition, the uncertainties caused 
by changing the signal function, background shape, and 
the fitting range in the fit of the invariant mass spectra 
recoil 7r"'"7r~ are also taken into account. To estimate 
the uncertainty caused by a change of the signal func- 
tion, we also fit the J /ip peak with the J /ip histogram 
shape, which is obtained from the recoil mass spectrum 
of vr+TT^ in ij:' — > 1:^71^ J/ ij), J /if) — >■ /i'^/i". The change 
of the result is just 0.04%. The uncertainty by changing 
the background shape from a second order Chebychev 
function to a first order one is less than 0.16%. If the fit- 
ting range is changed from [3.07, 3.13] GeV/c^ to [3.08, 
3.12] GeV/c^, the change is 0.32%. The total systematic 
uncertainty from the fitting, 0.37%, is the sum of these 
errors in quadrature. 



background normalized with and without the energy- 
dependent factor is taken as a background uncertainty. 

To estimate the background uncertainty from the 
beam-gas events, we select samples of beam-gas events 
in the J/?/' and continuum data. The candidate beam- 
gas events must have one or two charged tracks with the 
points of closest approach satisfying \Vz\ > 5 cm and 
\Vz\ < 15 cm and the visible energy less than 0.5 GeV. 
26844720 events are selected from the J/tp data, corre- 
sponding to 93470 events expected in the continuum data 
by normalizing with the integrated luminosities. Com- 
pared with 96230 beam-gas events directly selected from 
the continuum data, the difference between them, 3%, is 
taken as a background uncertainty. 

By adding all the above effects in quadrature, the total 
background uncertainty is 3.6%. Since the background 
ratio in J/t/j — )■ inclusive events is 3.5%, the systematic 
uncertainty in the number of J/tp events is 0.13%. 

E. Dependence on charged multiplicity 



D. Background uncertainty 

In the calculation of the number of J/^p events, the 
background events from QED processes, cosmic rays and 
beam-gas events are estimated by normalizing the se- 
lected continuum events by the integrated luminosities 
according to Eq. ([3]). Therefore the statistical error of 
the number of events selected from the continuum data, 
0.69% and the uncertainties due to the measurement of 
the integrated luminosities of the Jf-tjj data and contin- 
uum data, 1.5%, must be taken into account in the back- 
ground uncertainty. 

As discussed in Section 3, normalizing cosmic rays 
and beam-gas events with the energy-dependent factor 
for QED processes is not correct. To account for this, 
the difference, 1.1%, between the determinations of the 
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FIG. 8. The distributions of the charged multiplicity of 
J/tp inclusive events for J/ip data (dots with error bars) 
and tp' data (squares with error bars) and MC simulation 
generated with and without the Lundcharm model (solid and 
dashed histograms, respectively). 



In order to reduce the number of beam-gas events 
in this analysis, the selected J/ip — >■ inclusive events 
are required to have at least two good charged tracks 
{Ngood ^ 2). The uncertainty from this requirement is 
estimated by varying the charged multiplicity require- 
ment from Ngood > 2 to Ngood > 3. For comparison, the 
values obtained for the two cases are listed in Table |lTl 
The change of the number of J /ip events, 0.76%, is taken 
as the systematic uncertainty of the charged multiplicity 
requirement. 



F. Noise mixing 

Noise in the BESIII detector has been included in the 
realization of MC simulation by mixing in noise from 
events recorded using a random trigger for both J /ip and 
Ip' data. To determine the systematic error associated 
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FIG. 9. The distributions of deposited energy in EMC of 
photon in e^e~ 77 for data (dots) and MC simulation 
(histogram) . 
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with the noise reahzation in MC simulation, the tp' MC 
sample is reconstructed with the higher noise from J/tp 
data, and the change of the detection efficiency correction 
factor, 0.4%, is taken as a systematic uncertainty in the 
determination of the number of J/tp events. 

In this analysis 106 million of ip' events are used to 
determine the detection efhciency. However, the noise 
level was not entirely stable during the period of ip' data 
taking. To check the effect of the changing noise level on 
the detection efficiency, the ip' data and the MC sample 
are divided into three sub-samples, and the detection ef- 
ficiency is determined for each of the three samples. The 
change of the detection efficiency and the correction fac- 
tor lead to a change in the number of J/ip events. The 
maximum change, 0.28%, is taken as the systematic un- 
certainty associated with the changing noise levels. The 
total systematic uncertainty from the noise mixing effect 
is estimated to be 0.49% by adding the individual error 
contributions in quadrature. 

G. Estimation of Nj/^ with the sideband of 

The reliability of the determination of the number of 
J/ip events obtained from the above method is checked 
by applying another method entailing two different pro- 
cedures. One difference concerns the selection of inclusive 
events, which is essentially the same as in Section 2. ex- 
cept for the requirement on the track vertex position Vz 
along the beam direction. Here we determine the average 
position Vz of the charged tracks. The signal region for 
inclusive events is defined by \Vz\ < 4 cm. This require- 
ment is also applied in the determination of the detection 
efficiency and the correction factor. The Vz distribution 
is shown in Fig. 1141 

The second difference is in the background estimation. 
The numbers of background events from cosmic rays and 
beam-gas interactions are estimated from the Vz side- 
band, defined by 6 < iVzl < 10 cm. The subtraction of 
the sideband events from the events in the signal region 
removes the cosmic ray and beam-gas events. 



The sideband subtraction does not account for the 
QED background events since the Vz distribution is simi- 
lar to that of inclusive events from J/tp decays. However 
the continuum data allows us to estimate the contribu- 
tion of QED processes in the inclusive events selected 
from J/^p data. The same event selection is applied to 
the continuum data to select the QED events. After sub- 
tracting the cosmic rays and beam-gas events estimated 
with the same sideband method as for J/tp data, the 
amount of background events from the QED processes in 
the selected inclusive events is estimated by normalizing 
according to the integrated luminosities of the continuum 
and J/tp data according to Eq. ([3]). 

The same procedures have been used to determine the 
detection efficiency from ip' data and the correction fac- 
tor with MC samples. At last, the number of J/tp events 
is determined to 224.9 million. The change in the number 
of J/tp events with respect to the previous method dis- 
cussed in chapter 6 is 0.20% and is taken as a systematic 
uncertainty. 

H. Selection efficiency uncertainty of two soft pions 

According to a MC study, the selection efficiency of soft 
pions, Ett+tt-, recoiling against J/tp in tp' — >■ n^ir^J/tp 
depends on the multiplicity of the J/tp decay. To study 
its effect on the determination of the number of J/tp 
events, tp' -J> 7r+7r~(7r°7r°) J/V', J/tp -J> /.t+^~, 2(7r+7r~) 
events are selected from data and inclusive MC samples, 
and then re- weighting factors are determined for J/ tp de- 
caying into different multiplicities by comparing the cor- 
responding selection efficiency of soft pions between data 
and MC. The difference between the results with and 
without re-weighting, 0.34%, is taken as the uncertainty 
due to the selection efficiency uncertainty of the soft pi- 
ons in tp' TT^TT^ J/tp. 

The systematic uncertainties from different sources 
studied above are listed in Table Hill The total systematic 
uncertainty, 1.24%, is the sum of them added in quadra- 




cos 

FIG. 10. The cos 9 distributions of photon in e"'"e~ — > 77 for 
data (dots) and MC simulation (histogram). 
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FIG. 11. The TT momentum distributions from tp' data (dots 
with error bars) and MC simulation of tp' — >■ tt'^tt~ J/tp, 
J/tp — >■ n'^ n~ (histogram). 
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ture. 



VIII. SUMMARY 

Using J/ip — >■ inclusive events, the number of J/V' 
events collected with the BESIII detector in 2009 is de- 
termined to be 

iVj/^ = (225.3 ±2.8) X 106, (5) 

where the error is the systematic error and the statistical 
one is negligible. 
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FIG. 12. The invariant mass recoiling against selected vr tt^ 
pairs for 7/;' data. A clear peak corresponding to i])' — >■ 
■K^-K~ J/^j), J I'll) inclusive is seen. The curves are the re- 
sults of the fit described in the text. 
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FIG. 13. The invariant mass recoiling against selected ■K'^ir~ 
pairs for ip' data. Here, in addition to selection criteria on the 
pion pairs, the remaining portion of the event must satisfy the 
selection criteria for J/tp — >■ inclusive events. The curves are 
the results of the fit described in text. 



TABLE I. The values of different parameters used in the cal- 
culation and the resulting number of J/tp events. 

item value 

Nsei (174.28 ±0.01) X 10" 

Nbg (5.96 ± 0.04) X lO'' 

etrig 1-00 

0.7391 ±0.0002 

etc 0.7587 ± 0.0006 

e^it 0.7669 ± 0.0004 

fcor 1.0108 ±0.0009 

iVj/V, (225.30 ± 0.02) x 10" 



TABLE II. The number of J/xp events and values used in the 
calculation for Ngood > 2 and Ngood > 3. 

item Ngood > 2 Ngood > 3 

Nsei 174.28x10" 119.89x10" 

Nbg 5.96 x 10® 1.70 X lO'' 

etrig 1.00 1.00 

4'afa 0.7391 0.5050 

0.7587 0.5451 

ei(^ 0.7669 0.5620 

foor 1.0108 1.0310 

Nj/i, 



TABLE III. Summary of systematic uncertainties on the num- 
ber of J/tp events. 



Sources 


Relative error (%) 


MC model uncertainty 


0.49 


Tracking efficiency 


0.40 


Fitting of J/tp peak 


0.37 


Background uncertainty 


0.13 


Multiplicity requirement 


0.76 


Noise mixing 


0.49 


Sideband method 


0.20 


e^+^~ uncertainty 


0.34 


Total 


1.24 



-10 10 

Vjcm) 



FIG. 14. The distributions of the average z vertex of charged 
tracks, t4, for J/i/> data (dots with error bars), MC simulation 
of J 1^ inclusive plus continuum data (sohd histogram) 
and continuum data (dashed histogram). 



